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Observation of Xcj decays 
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Decays of Xc0,2 — > wo; are observed for the first time using a sample of 14.0 x 10® '0(25') events collected with 
the BESII detector. The branching ratios are determined to be B(xcO ^ i^i^) = (2.29 ± 0.58 ± 0.41) x 10"^ and 
B{xc2 -* ijJi^) = (1.77 ± 0.47 ± 0.36) x 10"^, where the first errors are statistical and the second systematic. The 



1. Introduction 

Exclusive quarkoninum decays provide an im- 
portant laboratory for investigating pcrturbative 
quantum chromodynamics. Compared with J/ij) 
and ip{2S) decays, one has much less knowledge 
on PC = + + XcJ decays. While a few exclu- 
sive decays of Xcj have been measured, many de- 
cay modes remain unknown. Current theoretical 
analyses of XcJ decays provide only a rough treat- 
ment of the color-octet wave function. For XcJ — * 
vector vector mode, so far only measurements of 
XcJ # [1] and XcJ ^ X*(892)0if* (892)0 p] 
are available with low statistics. Precise measure- 
ments for more channels will help in better un- 
derstanding the various mechanism [3,4] of XcJ 
decays and the nature of ^Pj cc bound states. 

Further, the decays of XcJ, especially XcO and 
Xc2, provide a direct window on glueball dy- 
namics in the 0"*""*" and 2++ channels since the 
hadronic decays may proceed via cc — > 55 — > 
qqqq. 

Recently, the branching ratio for XcO ^ 
/o(980)/o(980) [5] has been measured by the BES 
collaboration. In the present analysis, a search for 
XcO, 2 decaying into 7r+7r~7r07r+7r~7r° final states 
is carried out using 14 million ip{2S) events [6] 
accumulated at the upgraded BES detector (BE- 
SII). Signals of XcO and Xc2 decaying to to pairs in 
il>{2S) radiative decays are observed for the first 
time. 

2. The BES detector 

The Beijing Spectrometer (BES) is a conven- 
tional solcnoidal magnet detector that is de- 
scribed in detail in Ref . [7] ; BESII is the upgraded 
version of the BES detector [8]. A 12-layer vertex 
chamber (VC) surrounding the beam pipe pro- 
vides trigger and position information. A forty- 
layer main drift chamber (MDC), located radi- 
ally outside the VC; provides trajectory and en- 
ergy loss {dE/dx) information for charged tracks 
over 85% of the total solid angle. The m omen- 
tum resolution is (Tp/p = 0.017-\/l + {p in 
GeV/c), and the dE/dx resolution for hadron 
tracks is ~ 8%. An array of 48 scintillation coun- 
ters surrounding the MDC measures the time-of- 



flight (TOF) of charged tracks with a resolution of 
~ 200 ps for hadrons. Outside of the TOF coun- 
ters is a 12-radiation-length barrel shower counter 
(BSC) comprised of gas proportional tubes inter- 
leaved with lead sheets. This measures the en- 
ergies of electrons and photons over ~ 80% of 
the total solid angle with an energy resolution of 
aE/E = 22%/^/^ {E in GeV). Outside of the 
solenoidal coil, which provides a 0.4 Tesla mag- 
netic field over the tracking volume, is an iron flux 
return that is instrumented with three double lay- 
ers of counters that identify muons of momentum 
greater than 0.5 GeV/c. 

A GEANT3 based Monte Carlo (MC) program 
with detailed consideration of the detector perfor- 
mance (such as dead electronic channels) is used 
to simulate the BESII detector. The consistency 
between data and Monte Carlo has been care- 
fully checked in many high purity physics chan- 
nels, and the agreement is quite reasonable [9] . 

3. Event selection 

3.1. WO) signal 

In this analysis, XcJ '^'^ 7r+7r~7r''7r+7r~7r° 
channels are investigated using ipi^S) radiative 
decays to XcJ- Events with four charged tracks 
and five or six photons are selected. Each 
charged track is required to be well fit by a three- 
dimensional helix and to have a polar angle, 9, 
within the fiducial region | cos 9\ < 0.8. To ensure 
tracks originate from the interaction region, we 

require Ky = ^jyf+y^ < 2 cm and \Vz\ < 20 
cm, where Vx, Vy, and Vz are the x,y and z co- 
ordinates of the point of closest approach of each 
track to the beam axis. 

A neutral cluster is considered to be a photon 
candidate if it is located within the BSC fiducial 
region (| cosd\ < 0.8), the energy deposited in the 
BSC is greater than 40 MeV, the first hit appears 
in the first 10 radiation lengths, and the angle 
between the cluster and the nearest charged track 
is greater than 6°. 

A six constraint (6-C) kinematic fit to the hy- 
pothesis ^j{2S) 77r^~7r~7r07r"'"7r~7rO with the in- 
variant mass of the two photon pairs constrained 
to the TT^ mass is performed, and the of the 6-C 
fit is required to be less than 15. For events with 
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six photons candidates, the combination having 
the minimum is chosen, and the probability of 

the 6-C fit is required to be larger than that of the 
7-C fit to the hypothesis V(25') 27r+27r-37r° 
to suppress potential background from tp{2S) 

Since there are four lj pair combinations from 
7r+7r~7r°7r"'"7r~7r°, the u) pair with the minimum 
R, which is defined as 

R = - 0-^83)2 + (M^\,-,„ - 0.783)2, 

is chosen for further analysis. Here, M^+^-^o is 
the invariant mass of three pions and superscript 
1, 2 denote different pion combinations. There- 
fore, there is only one entry for each event. 

Figures 1 and 2 show mass distributions for 
candidate events in the high mass {Mg^ > 3.2 
GeV/c^) and low mass regions {Mq„ < 3.2 
GeV/c^), respectively. Here (a) is the scatter 
plot of M^+^-^o versus M^+^-^o, (b) is the 
M^+^-^o distribution recoiling against the oppo- 
site LO, selected by requiring |M^+^-^o — 783| < 50 
MeV/c^, and (c) is the M^iu, invariant mass dis- 
tribution for events in the w pair signal region, 
defined by /2 < 50 MeV/c^. In Fig. 1, clear u! Sig- 
nal can be seen in (b), and clear Xco £md Xc2 sig- 
nals in (c), indicating the existence of Xco.2 — > ^^^^ 
decays. By contrast, in the low Mg^ mass region, 
shown in Fig. 2, the u) pair signal is less signifi- 
cant than in the high mass region. Here, only ui 
pair events in high mass region are studied. 

In order to test if the selection criteria in this 
analysis will give 'fake' u pair events from non-o; 
pair events, 300000 MC simulated tp IXco 
jGir events are generated in which XcO Stt de- 
cays according to the phase space. Fig. 3 shows 
the Mj^+^-^o distributions of the surviving MC 
phase space events after requiring the same selec- 
tion criteria as for the real data. No peak around 
the ui mass is seen, which shows that the iv pair 
selection criteria in this analysis does not gener- 
ate fake Lu pair signals. 

The annular region around the lo pair signal 
circle, shown in Fig 1(a), is taken as the side- 
band region. Fig. 4 shows the Me^ sideband dis- 
tributions defined using the radius R to be (a) 
150 < i? < 300 MeV/c^ and (b) 100 < i? < 200 
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Figure 1. Distributions of events surviving the se- 
lection criteria described in the text with MgTr > 



(a) 



versus 



3.2 GeV/c2. 
(b) M^^+T^-^o recoiling against the opposite u, 
selected by requiring |M^+„-„o — 783| < 50 
MeV/c^, and (c) M,^^ invariant mass distribu- 
tion for events where the w pair satisfies ii < 50 
MeV/c2. 



MeV/c^. No obvious XcJ signals seen in these 
sideband distributions. 

3.2. MC simulation 

A MC simulation of ip{2S) 7XcJ) XcJ '^'^ 
is used to determine the detection efficiency. The 
proper angular distributions of the photon emit- 
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phase space simulated iIj{2S) JXco, Xco 
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Figure 2. Distributions defined as in Fig. 1 but 
with Mevr < 3.2 GeV/c^ 



Figure 4. MgTr distribution of events in sideband 
regions (a) 150 < i? < 300 MeV/c^ and (b) 100 < 
ii< 200 MeV/c^. 



ted in ipCiS) — > 7XcJ are used [10]. Fig. 5 shows 
the distributions, identical to those in Fig. 1 for 
MC simulated ip{2S) — » 7Xco, Xco events 
passing the same selection criteria as for the real 

data. MC simulated '0(25') 7Xc2, Xc2 — > i^i^ 
events have similar distributions. 

3.3. Mass spectrum fit 

The Maximum Likelihood (ML) method is used 
to fit the M^^jj mass spectrum of events in the 
U! pair signal region (Fig. 1(c)). The xo,2 signal 
functions are determined from MC simulation, as 
shown in Fig. 5(c) for XcO) while the background 



function is taken from the sideband distribution, 
shown in Fig 4(a). The fit result is represented 
by the solid curve in Fig. 6, and the fit yields 

AT^^o = 38.1 ± 9.6, 7Vx=2 = 27.7 ± 7.4. 

The statistical significances of XcO and Xc2 are 
4.4(7 and 4.7(t, respectively, which arc estimated 
from v2Aln£, where Aln£ is the difference be- 
tween the logarithmic ML values of the fit with 
and without the corresponding signal function. 
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Figure 5. Distributions defined as in Fig. 1 from 
MC simulated ipi'^S) "JXcO, XcO events. 



4. Systematic error 

Tlie systematic error in tliis branching ratio 
measurement includes the uncertainties in the 
MDC tracking efficiency, photon efficiency, kine- 
matic fit, background shape, number of '0(25') 
events, etc. 

4.1. MDC tracking efficiency and photon 
efficiency 

For charged tracks, the uncertainty of the 
tracking efficiency is determined by comparing 
data and MC for some very clean J/ i/j decay chan- 
nels [9] , and an error of 2% is found for each track. 




(GeV/c ) 



3.70 



Figure 6. Fit of the M^^ distribution. Dots with 
error bars are data, the solid histogram represents 
the maximum likelihood fit result, and the dashed 
histogram is the sideband background. 



A similar comparison has also been performed for 
photons [12], and the difference is also about 2% 
for a single photon. 

4.2. Kinematic fit 

The systematic error associated with the kine- 
matic fit is due to differences between data and 
MC simulation in the determination of the track 
momentum, the track fitting error matrix, and 
the photon energy and direction. The effect is 
studied for charged tracks and neutral tracks sep- 
arately. By comparing the number of events be- 
fore and after the kinematic fit for very clean 
event samples for data and MC simulated data, 
the difference is determined to be 8.4%, which is 
taken as the systematic error. 

4.3. Background shape 

Two different sideband Mqt^ spectrum shapes, 
shown in Fig. 4, are used as the background func- 
tion. The difference in the number of Xc0,2 events 
obtained with the two different shapes is taken as 
a systematic error. 

4.4. Binning, fit range, and signal region 

The differences caused by different binning and 
fit ranges in the ujlj mass spectrum fit are 1.2% 
and 3.4% for XcO and Xc2, respectively. Differ- 
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Table 1 

Individual sources and total systematic error (%). 



Source 


XcO tJW 


Xc2 U)U 


track efficiency 


8 


8 


photon efficiency 


10 


10 


6-C fit 


8.4 


8.4 


background shape 


6.0 


1.0 


signal region 


3.4 


4.3 


binning and fit range 


1.4 


3.2 


angular distribution 




9.4 


No. of ip{2S) events 


4 


4 


B{^{2S) ^ 7Xcj) 


5.1 


6.7 


B{uj 3tt) 


0.9 


0.9 


B(7rO ^ 77) 


0.0 


0.0 


Total 


18.1 


20.4 



ent sized signal regions yield differences of 3.1% 
and 2.4% for XcO and Xc2, respectively, which are 
taken as a systematic error. 



5. Results 

The branching ratio of B{XcJ 
mined from 



u)u) is deter- 



bs 



i,{2S) 



h 



fi 



fi 



where N°''' 

XcJ 



is the number of events selected, 
N^(2S) the total number of tp{2S) events, e is 
the detection efficiency for the investigated chan- 
nel, and /i , /2 and /a are the branching ratios of 
ip{2S) jXcJ,^^ — > Stt, and 7r° — » 77, respec- 
tively. Table 2 lists the XcO,2 branching 
ratio results, together with numbers used in the 
branching ratio calculation. 

In summary, lolo signals in the decay of Xc0,2 
are observed, and their branching ratios mea- 
sured for the first time. XcO and Xc2 decays to 
LOLU have similar decay branching ratios, which is 
different from other Xcj VV decays {XcJ — * 
00, i^*(892)0K*(892)°). This measurement, to- 
gether with previous measurements of Xcj — * 
VV, will be helpful in understanding the nature 
of XcJ states. 



4.5. Angular distribution of Xcj — >■ wo; 

In the estimation of the efficiency, a phase space 
generator with only the angular distribution of 
the radiative photon is considered. While this is 
correct for XcO decays, it may introduce bias for 
Xc2 decays. The effect is estimated by generat- 
ing different angular distributions of the omega 
in the Xc2 rest frame. The efficiency difference 
between these tests and the phase space genera- 
tor is estimated to be 9.4%, which is put into the 
systematic error. 

4.6. Branching ratios of intermediate 
states 

The errors on intermediate state branching 
ratios are obtained from the PDG [13] except 
for B{i>{2S) IXcj), where recent CLEO re- 
sults [14] arc used. Table 1 summarizes all con- 
tributions to the systematic errors, and the total 
systematic error is determined by the quadratic 
sum of all terms. 
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